New results for the two neutrino double beta decay in deformed 
nuclei with angular momentum projected basis 

A. A. Raduta^' 6 ), A. Escuderos^ , Amand Faessler d \ E. Moya de Guerra c - ) and P.Sarriguren c ) 
Institute of Physics and Nuclear Engineering, Bucharest, POBox MG6, Romania 
^Department of Theoretical Physics and Mathematics, 
Bucharest University, POBox MG11, Romania 
^Institute de Estructura de la Materia, 



O 

o . 

fSj . Consejo Superior de Investigationes Cientificas, 



Serrano 123, E-28006 Madrid, Spain and 

O ' Institut fuer Theoretische Physik, Universitaet Tuebingen, D-72076 Tuebingen, Germany 

(N 

Abstract 

00 1 Four nuclei which are proved to be 2z//?/3 emitters ( 76 Ge, 82 Se, 150 Nd, 238 U), and four suspected, 

(N ■ 

due to the corresponding Q-values, to have this property ( 148 Nd, 154 Sm, 160 Gd, 232 Th), were 

t-h ; 

treated within a proton-neutron quasiparticle random phase approximation (pnQRPA) with a 

^ ■ 

projected spherical single particle basis. The advantage of the present procedure over the ones using 
a deformed Woods Saxon or Nilsson single particle basis is that the actual pnQRPA states have a 
definite angular momentum while all the others provide states having only K as a good quantum 
number. The model Hamiltonian involves a mean field term yielding the projected single particle 
states, a pairing interaction for alike nucleons and a dipole-dipole proton-neutron interaction in 
both the particle-hole (ph) and particle-particle (pp) channels. The effect of nuclear deformation 
on the single beta strength distribution as well as on the double beta Gamow- Teller transition 
amplitude (Mqt) is analyzed. The results are compared with the existent data and with the 
results from a different approach, in terms of the process half life Ti/2- The case of different 
deformations for mother and daughter nuclei is also presented. 

PACS numbers: 23.40.Hc, 21.60.Jz, 27.50. +c 
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I. INTRODUCTION 



One of the most exciting nuclear physics subject is that of double beta decay. The 
interest is generated by the fact that in order to describe quantitatively the decay rate one 
has to treat consistently the neutrino properties as well as the nuclear structure features. 
The process may take place in two distinct ways: a) by a 2vj3j3 decay the initial nuclear 
system, the mother nucleus, is transformed in the final stable nuclear system, usually called 
the daughter nucleus, two electrons and two anti-neutrinos b) by the Qvf}(3 process the final 
state does not involve any neutrino. The latter decay mode is especially interesting since 
one hopes that its discovery might provide a definite answer to the question whether the 
neutrino is a Majorana or a Dirac particle. The Qvj3j3 decay is an extremely rare process and 
moreover it is hard to distinguish the electrons emerging from the two processes. For some 
processes there exists information about the low limits of the process half-lives. Combining 
this information with the nuclear matrix elements, some conclusions about the upper limits 
of both neutrino effective mass and effective right-handedness of the electroweak interaction 
was possible. Unfortunately there are no reliable tests for the nuclear matrix elements 
involved and therefore some indirect methods should be adopted. It is worth mentioning 
that the matrix elements which are responsible for neutrinoless double beta decay are similar 
to those needed for calculating the 2vj3f3 decay rate, for which there exists experimental data. 
Due to this feature an indirect test for the matrix elements used for 0^/3/3 is to use those 
m.e. which describe quantitatively the 2v(3(3 decay. 

For such reasons many theoreticians focused their efforts in describing consistently the 
data for 2vf3f3 decay. The contributions over several decades have been reviewed by many 
authors. Instead of enumerating the main steps achieved toward improving the theoretical 
description we advise the reader to consult few of the review works |2j, mag. 

It is interesting to note that although none of the double beta emitters is a spherical 
nucleus most formalisms use a single particle spherical basis. More than 10 years ago, two 
of us |9] proposed a formalism to describe the process of two neutrinos double beta decay 
in a projected spherical basis. It was for the first time that a pnQRPA approach for a two 
body interaction in the ph and pp channels with a deformed single particle basis was per- 
formed. Moreover, effects which are beyond the proton-neutron quasiparticle random phase 
approximations (pnQRPA) have been accounted for by means of a boson expansion proce- 
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dure. A few years later the influence of nuclear deformation upon the contribution of the 
spin-flip configurations to the Gamow- Teller double beta transition amplitude, was studied 



10j . In the meantime several papers have been devoted to the extension of the pnQRPA 



procedure to deformed nuclei, the applications being performed for studying the single beta 
decay properties as well as the double beta decay rates. Thus, pnQRPA approaches us- 
ing as deformed single particle basis, Nilsson or deformed Woods Saxon states have been 



formulated 
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12l Il3j . Also a selfconsistent deformed method was formulated where the 
single particle basis was obtained as eigenstates of a deformed mean field obtained thr ough 
a Hartree-Fock treatment of a density dependent two body interaction of Skyrme type (12j ■ 
The present investigation is, in fact, a continuation of the work from Ref. [9j. Therein the 
single particle energies were depending linearly on a parameter which simulates the nuclear 
deformation. By contrast, here the core volume conservation constraint, ignored in the 
previous paper, determines a nonlinear deformation dependence for single particle energies. 
Of course, having different single particle energies one expects that the pairing properties 
and the double beta matrix elements are modified. Another issue addressed in the present 
paper is whether considering different deformations for the mother and daughter nuclei, 
modifies significantly the double beta transition amplitude (Mgt). To be more specific, 
we recall that the standard pnQRPA approach including only the two body interaction in 
the particle-hole (ph) channel yields a Mqt value much larger than the experimental value 
extracted from the corresponding half-life. Apparently, the desired Mqt suppression might 
be obtained by a suitable choice of the two body interaction in the particle-particle (pp) 
channel. However, the fitted strength is close to the value where Mgt cancels and moreover 
close to the critical value where the pnQRPA breaks down. It is obvious that increasing 
the deformation for the daughter nucleus the pnQRPA phonon state is less correlated and 
therefore the pnQRPA breaking point is pushed toward larger values. In this respect, one 
may say that the value of the pp interaction strength which reproduces the experimental 
value for Mgt becomes reliable, i.e. the corresponding pnQRPA ground state of the daughter 
nucleus is stable against adding anharmonic effects. 

The formalism and results of the present paper will be presented according to the fol- 
lowing plan. In Section II a brief review of the projected spherical single particle basis will 
be presented. Section III deals with the pnQRPA treatment of a many body Hamiltonian 
which describes the nuclear states of the mother, daughter and intermediate odd-odd nuclei, 
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involved in the 2v(3(3 process. In Section IV, we discuss the results for eight double beta 
emitters: 76 Ge, 82 Se, 148 Nd 150 Nd, 154 Sm, 160 Gd, 232 Th, 238 U for which the strength distri- 
bution for single (3~ and f3 + emission for mother and daughter nuclei, respectively, the Mqt 
and half lives values for the double beta decay process are presented. A short summary and 
concluding remarks are given in Section V. 



II. PROJECTED SINGLE PARTICLE BASIS 

In Ref. one of us, (A.A.R.), introduced an angular momentum projected single 

particle basis which seems to be appropriate for the description of the single particle motion 
in a deformed mean field generated by the particle-core interaction. This single particle 
basis has been used to study the collective Ml states in deformed nuclei 3] as well as the 
rate of double beta process 0, [h]] . Recently a new version has been proposed where the 
deformation dependence of single particle energies is nonlinear and therefore more realistic 
jlfil Ir^ . In order to fix the necessary notations and to be self-contained, in the present 
work we describe briefly the main ideas underlying the construction of the projected single 
particle basis. Also some new properties for the projected basis are pointed out. 

The single particle mean field is determined by a particle-core Hamiltonian: 

H = H sm + H core -Mu 2 r 2 £ Yl «1/a„. (2.1) 

A=0,2 -\<ti<\ 

where H sm denotes the spherical shell model Hamiltonian while H core is a harmonic 
quadrupole boson (6+) Hamiltonian associated to a phenomenological core. The interac- 
tion of the two subsystems is accounted for by the third term of the above equation, written 
in terms of the shape coordinates a 00 , a 2fl - The quadrupole shape coordinates are related to 
the quadrupole boson operators by the canonical transformation: 

a 2 » = ^(bl+(-rb 2 ^), (2.2) 

where k is an arbitrary C number. The monopole shape coordinate is to be determined from 
the volume conservation condition. In the quantized form, the result is: 

1 



5+E(2^ + (^ t + & -A)(-r) 



(2.3) 
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Averaging H on the eigenstates of H sm , hereafter denoted by \nljm), one obtains a deformed 
boson Hamiltonian whose ground state is, in the harmonic limit, described by a coherent 
state 

y g = exp[d(b+ -b 20 )}\0) b , (2.4) 

with |0)& standing for the vacuum state of the boson operators and d a real parameter which 
simulates the nuclear deformation. On the other hand, the average of H on is similar 
to the Nilsson Hamiltonian Due to these properties, it is expected that the best trial 
functions to generate, through projection, a spherical basis are: 

= \nljm)* g . (2.5) 

The upper index appearing in the l.h. side of the above equation suggests that the product 
function is associated to the particle-core system. The projected states are obtained, in the 
usual manner, by acting on these deformed states with the projection operator 

P mk = J & MK m (fl)R(fl)dSl. (2.6) 

We consider the subset of projected states : 

*3j (d) = KiiPUlnljl}^] = (d). (2.7) 

which are orthonormalized. 

The main properties of these projected spherical states are: a) They are orthogonal with 
respect to I and M quantum numbers, b) Although the projected states are associated 
to the particle-core system, they can be used as a single particle basis. Indeed, when a 
matrix element of a particle like operator is calculated, the integration on the core collective 
coordinates is performed first, which results in obtaining a final factorized expression: one 
factor carries the dependence on deformation and one is a spherical shell model matrix 
element, c) The connection between the nuclear deformation and the parameter d entering 
the definition of the coherent state (2.4) is readily obtained by requiring that the strength 
of the particle-core quadrupole-quadrupole interaction be identical to the Nilsson deformed 
term of the mean field: 
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Here Q± and Q z denote the frequencies of Nilsson's mean field related to the deformation 



5 = ^45/16^/3 by: 



n± = (^)^, n z = (f^r 2/3 . (2.9) 



The constant k was already defined by Eg 12. 21 This is at our disposal since the canonical 
property of the quoted transformation is satisfied for any value of k. The average of the 
particle-core Hamiltonian H' = H — H core on the projected spherical states defined by 
Eq.(2.7) has the expression 



+ fuv (N+- 



1 + ^7 + 
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Here we used the Condon-Shortley convention and notation for the Clebsch Gordan coeffi- 
cients m . Ij stands for the following integral 

if = f 1 Pj(x)[P 2 (x)} i exp\d 2 P 2 (x)]dx, i = 0, 1. (2.11) 
Jo 

where Pj(x) denotes the Legendre polynomial of rank J. It is worth mentioning that the 
norms for the core's projected states as well as the matrix elements of any boson operator on 
these projected states can be fully determined once the overlap integrals defined in (|2.11jl . 
are known ^|. Since the core contribution does not depend on the quantum numbers of 
the single particle energy level, it produces a shift for all energies and therefore is omitted in 
Eq.(2.10). However, when the ground state energy variation against deformation is studied, 
this term must be included. 

The first term from ()2.10j) is, of course, the single particle energy for the spherical shell 
model state \nljm). The second term, linear in the deformation parameter d, is the only one 
considered in the previous works devoted to the double beta decay of deformed nuclei within 
a projected spherical basis formalism. The third term from Eq. (j2.10j) is determined by the 
monopole-monopole particle-core coupling term after implementing the volume conservation 
condition. This term is the one responsible for the nonlinear deformation dependence of e^,-. 
The energies are represented as function of the deformation parameter d, for the major 
shells with N equal to 3 and 4, in Fig. 1. We remark that the energies shown in the above 
mentioned plot depend on deformation in a different manner than those obtained in Ref. 
Q|. Indeed, therein they depend linearly on deformation, while here non-linear effects are 



present. The difference between the two sets of energies is caused by the fact that here 
the volume conservation condition was used for the monopole shape coordinate, while in 
Ref.[l^| this term is ignored. The difference in the single particle energies is expected to 
cause significant effects on the single and double beta transition probabilities. Actually, this 
is the main motivation for the present investigation. 

As shown in Fig. 1, the dependence of the new single particle energies on deformation is 
similar to that shown by the Nilsson modelQ]. 

Although the energy levels are similar to those of the Nilsson model, the quantum numbers 
in the two schemes are different. Indeed, here we generate from each j a multiplet of (2j + 1) 
states distinguished by the quantum number I, which plays the role of the Nilsson quantum 
number Q, runs from 1/2 to j and moreover the energies corresponding to the quantum 
numbers K and -K are equal to each other. On the other hand, for a given I there are 21 + 1 
degenerate sub-states while the Nilsson states are only double degenerate. As explained in 
Ref . [yj , the redundancy problem can be solved by changing the normalization of the model 
functions: 

(<^ M lO = i =}► £($i A/ |^ M > = 2 - (2-12) 

M 

Due to this weighting factor the particle density function is providing the consistency result 
that the number of particles which can be distributed on the (21+1) sub-states is at most 
2, which agrees with the Nilsson model. Here a stands for the set of shell model quantum 
numbers nlj. Due to this normalization, the states used to calculate the matrix elements 



of a given operator should be multiplied with the weighting factor J 2/(21 + 1). The role 
of the core component is to induce a quadrupole deformation for the matrix elements of 
the operators acting on particle degrees of freedom. Indeed, for any such an operator the 
following factorization holds: 

(^•II^H^W) = fSrinimiln'l'f). (2.13) 

The factor / carries the dependence on the deformation parameter d while the other factor 
is just the reduced matrix elements corresponding to the spherical shell model states. For 



details we advise the reader to consult Refs. }l4 



Q. 



Concluding, the projected single particle basis is defined by Eq. (|2.7|) . Although these 
states are associated to a particle-core system, they can be used as a single particle basis 
due to the properties mentioned above. 



Therefore, the projected states might be thought of as eigenstates of an effective rotational 
invariant fermionic one-body Hamiltonian H e ff, with the corresponding energies given by 

Eq. (jrruD . 

He ff K M = tlml M . (2.14) 

This definition should be supplemented by the request that the matrix elements of any 
operator between states and $„/ M ', are given by Eq. (j2.13|) . Due to these features, 
these states can be used as single particle basis to treat many body Hamiltonians which 
involve one-body operators. This is the case of Hamiltonians with two body separable 
forces. As a matter of fact, such a type of Hamiltonian is used in the present paper. 

According to our remark concerning the use of the projected spherical states for describing 
the single particle motion, the average values may be viewed as approximate expressions 
for the single particle energies in deformed Nilsson orbits We may account for the 
deviations from the exact eigenvalues by considering, later on, the exact matrix elements of 
the two body interaction when a specific treatment of the many body system is applied. 

Few words about the vibrational limit, d — > 0, for the projected basis are necessary. It 
can be proved that the following relations hold: 

gm(A^) _1 = S Idt 
lim(^.||T fe ||$^,) = {nlj\\T k \\n'l'j'), 

l }™ e ilj = e nlj- (2-15) 

Note that in the limit d — *■ 0, the norms of the states with I ^ j are not defined while the 
limit of the I=j state, normalized to unity, is just the product state \nljM)\0) b . Indeed, 
the fourth equation ()2.15|) is fulfilled by neglecting a small quantity (g^p\) caused by the 
zero point motion term of the monopole-monopole particle-core interaction. Although in the 
limit d — > the norm of the states / ^ j is not defined, the limit of with I ^ j, exists. 
However, the corresponding energies are not identical to but very close to the spherical shell 
model state energy. 

lim ^ = + h M N + §) (§ + i (j - I + 1(1 - (-)-'))) Jp, j * I. (2.16) 
Indeed, this term should be compared with the spherical oscillator energy ({N + ^)Huq) from 
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e n ij. Since the factor 



i 



is very small (see Table 1) the correction of the shell model term 



l-A' 2 



is negligible. 

Due to the properties mentioned above, we may state that in the vibrational limit, d — > 0, 
the projected spherical basis goes to the spherical shell model basis. 

To complete our description of the projected single particle basis, we recall a fundamental 
result obtained in Ref.[l7|, concerning the product of two single particle states which com- 
prises a product of two core components. Therein we have proved that the matrix elements 
of a two body interaction corresponding to the present scheme are very close to the matrix 
elements corresponding to spherical states projected from a deformed product state with 
one factor as a product of two spherical single particle states, and a second factor consisting 
of a common collective core wave function. The small discrepancies of the two types of 
matrix elements could be washed out by using slightly different strengths for the two body 
interaction in the two methods. 

III. THE MODEL HAMILTONIAN AND ITS PNQRPA APPROACH 

As we already stated, in the present work we are interested to describe the Gamow- Teller 
two neutrino double beta decay of an even-even deformed nucleus. In our treatment the 
Fermi transitions, contributing about 20% and the "forbidden" transitions are ignored which 
is a reasonable approximation for the two neutrino double beta decay in medium and heavy 
nuclei. Customarily, the 2vf3f3 process is conceived as two successive single (3~ transitions. 
The first transition connects the ground state of the mother nucleus to a magnetic dipole 
state 1 + of the intermediate odd-odd nucleus which subsequently decays to the ground state 
of the daughter nucleus. Going beyond the pnQRPA procedure by means of the boson 
expansion procedure we were able to consider the process leaving the final nucleus in an 
excited collective state [19]. Such processes are not treated in the present paper. The states, 
mentioned above, involved in the 2v(3f3 process are described by the following many body 
Hamiltonian: 



H 




(3.1) 
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The operator c| q/m ( c ™/a/) creates (annihilates) a particle of type r (=p,n) in the state 
$^ M , when acting on the vacuum state |0). In order to simplify the notations, hereafter 
the set of quantum numbers a(= nlj) will be omitted. The two body interaction consists 
of three terms, the pairing, the dipole-dipole particle hole (ph) and the particle-particle 
(pp) interactions. The corresponding strengths are denoted by G T ,x,Xi, respectively. All 
of them are separable interactions, with the factors defined by the following expressions: 



2 



PtI — E o T , i C rIM c 



t J 



2T + 1 T1M tIM 7 

M * l ~ L 



Piiipn) = E ^(pIM\a^\nI'M')^3-cl IM c nrM >, 

MM' 1 1 



P{M) = E ^m*M'M')^cl IM cl~. (3.2) 

M,M' 1 1 

The remaining operators from Eq. (|3.1|) can be obtained from the above defined operators 
by hermitian conjugation. 

The one body term and the pairing interaction terms are treated first through the 
standard BCS formalism and consequently replaced by the quasiparticle one body term 
J2riM E T a\. IM a T iM- In terms of quasiparticle creation (a TlM ) and annihilation (cl t im) oper- 
ators, related to the particle operators by means of the Bogoliubov-Valatin transformation, 
the two body interaction terms, involved in the model Hamiltonian, can be expressed just 
by replacing the operators (3.2) by their quasiparticle images: 

p;(k) = <7 fc Ai M (A0 + **A w (A;)H^^ 

ft(k) = -[*fc4#) + ^ w (*x-) 1_M ^ , 

P+(k) = _[_^4 m (A;)+ % A^ (3.3) 

In the above equations the argument "k" stands for the proton-neutron state (p,n). Here, 
the usual notations for the dipole two quasiparticle and quasiparticle density operator have 
been used: 
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B\fj,(pn) — ^2 Cmp -m„ n a plpmp a nlnmn( Y" m "- (3-4) 
m p ,m n 

The coefficients a and 77 are simple expressions of the reduced matrix elements of the 
Pauli matrix a and U and V coefficients: 

2 2 

J-'n J-'n 

2 2 

^fc = ^C^pNI-U^A,,, = lyifplMfnWlpVln, ( 3 - 5 ) 

"ra "to 

The quasiparticle Hamiltonian is further treated within the proton-neutron random phase 
approximation (pnQRPA), i.e. one determines the operator 

Tl M = J2[X(k)Al(k) - Y(k)A x ^{k)(-f-% (3.6) 
k 

which satisfies the restrictions: 



Fim, rI M ,] = [# gp , rg = (3.7) 

These operator equations yield a set of algebraic equations for the X (usually called 
forward going) and Y (named back-going) amplitudes: 



{-B -a)( X y)A X y)' <38) 

Y\\X (*)| 2 - |i'(<;)| 2 ] = 1. (3.9) 

k 

The pnQRPA matrices A and B have analytical expressions: 



A k ,k' = {E p + E n )5 ppl 5 nn , + 2x(o k o k , + a k a k >) - 2xi{VkVk> + VkVk>), 

&k,k> = 2x{d k a k > + a k a k >) + 2xi{fj k r]k' + (3-10) 

All quantities involved in the pnQRPA matrices have been already defined. Note that the 
proton and neutron quasiparticle energies are denoted in an abbreviated manner by E p and 
E n , respectively. 

As can be seen from Eq. (|3.1j) the ph interaction is repulsive while the pp interaction 
has an attractive character. Due to this feature, for a critical value of xi the lowest root of 
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the pnQRPA equations may become imaginary. Suppose that xi is smaller than its critical 
value and therefore all RPA solutions (i.e. u) are real numbers and ordered as: 

U\ < ^2 < ••• < un s - (3-11) 

Here N s stands for the total number of the proton-neutron pair states whose angular mo- 
menta can couple to 1 + and moreover their quantum numbers n, I are the same. Hereafter 
the phonon amplitudes X and Y will be accompanied by a lower index "i" suggesting that 
they correspond to the energy oji. 

Since our single particle basis states depend on the deformation parameter d, so do the 
pnQRPA energies and amplitudes. The pnQRPA ground state (the vacuum state of the RPA 
phonon operator) describes an even-even system which might be alternatively the mother 
or the daughter nucleus. In the two cases the gauge and nuclear deformation properties are 
different which results in determining distinct pnQRPA phonon operators acting on different 
vacua describing the mother and daughter ground states, respectively. Therefore, one needs 
an additional index distinguishing the phonon operators of the mother and daughter nuclei. 
The single phonon states are defined by the equations: 

I';,''),-./ '•/: k = l,2,...N a . (3.12) 

Here the indices % and / stand for initial (mother) and final (daughter) nuclei, respectively. 
This equation defines two sets of non-orthogonal states describing the neighboring odd-odd 
nucleus. The states of the first set may be fed by a beta minus decay of the ground state 
of the mother nucleus while the states of the second set are populated with a beta plus 
transition operator from the ground state of the daughter nucleus. 

If the energy carried by leptons in the intermediate state is approximated by the sum of 
the rest energy of the emitted electron and half the Q- value of the double beta decay process 

AE = m e c 2 + ~Qpp, (3.13) 

the reciprocal value of the 2v/3(3 half life can be factorized as: 

(T 2^ r i = F|M GT (0+ - 0+)| 2 , (3.14) 

where F is an integral on the phase space, independent of the nuclear structure, while Mqt 
stands for the Gamow- Teller transition amplitude and has the expression j^: 
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m gt = VsJ2 



(oii^im>tt(i fc |iy>//(iyii)9;iio)/ 



(3.15) 



kk/ 



E k + AE + E 1+ 



In the above equation, the denominator consists of three terms: a) AE, which was already 
defined, b) the average value of the k-th pnQRPA energy normalized to the particular value 
corresponding to k=l, i.e. 



and c) the experimental energy for the lowest 1 + state. The indices carried by the transition 
operators indicate that they act in the space spanned by the pnQRPA states associated to 
the initial (i) or final (f) nucleus. Details about the overlap matrix of the single phonon 
states in the mother and daughter nuclei are given in Appendix A. 

Before closing this section we would like to say a few words about what is specific to our 
formalism. As we mentioned before the pnQRPA matrices depend on the deformation pa- 
rameter and therefore the RPA energies and states depend on deformation. Moreover, in the 
case that the mother and daughter nuclei are characterized by different nuclear deformations 
the RPA output for the two nuclei are affected differently by deformation. These features 
make the pnQRPA formalism build up with a deformed single particle basis quite tedious. 
Besides these difficulties, one should keep in mind the fact that the usual approaches define 
the states from the intermediate odd-odd nucleus not as a state of angular momentum 1 but 
states of a definite K, i.e., K = ±1,0. Under such circumstances from the pnQRPA states, 
the components of good angular momentum are to be projected out. This operation is usu- 
ally performed in an approximative way, (by transposing the result obtained in the intrinsic 
frame, to the laboratory frame of reference) which might be justified only in the strong 
coupling regime. Unfortunately, the double beta emitters are only moderately deformed, 
which makes the approximation validity, questionable. Actually this is the reason why the 
answer to the question of how much the results obtained with deformed single particle basis 
differ from the ones obtained with projected many body RPA state, is not yet known. 

By contrast, since our single particle states are projected spherical states, the RPA for- 
malism is fully identical to that one which is usually employed for spherical nuclei. Since in 
the vibrational limit, (d — > 0), our basis goes to the spherical shell model basis, one may say 
that the present formalism provides a unified description of spherical and deformed nuclei. 




(3.16) 
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IV. NUMERICAL RESULTS 



The formalism described in the previous sections was applied to eight nuclei among which 
four are proved to be double beta emitters ( 76 Ge, 82 Se, 150 Nd, 238 U) [22j and four suspected, 
due to the corresponding Qpp value, to have this property. 

The spherical shell model parameters are those given in Ref. i.e. 

huj = UA- 1/3 , C = -2huj K, D = -htoofi. (4.1) 

For the proton system, the pair of strength parameters (k, fi) takes the values (0.08;0.) for 
76 Ge, 76 Se, 82 Se, 82 Kr, (0.0637;0.6) for 148 > 150 Nd, 148 " 154 Sm, 154 ' 160 Gd, 160 Dy, (0.0577;0.65) for 
232 Th, 232 > 238 U, 238 Pu, while for the neutron systems of the three groups of nuclei mentioned 
above, the values are (0.08;0.), (0.0637,0.42), (0.0635;0.325), respectively. 

The projected spherical single particle basis, used in our calculations, depend on another 
two parameters, the deformation d and the factor k of the transformation (2.2) relating 
the boson operators with the quadrupole collective coordinate, according to Eq. (|2.2|) . 
These were fixed as follows. For the lightest nuclei, Ge, Se, Kr, involved in the process of 
the double beta decay, the two parameters were taken so that the relative energies of the 
states |1/||) and |lrf||) as well as the lowest root of the pnQRPA equations with a QQ 
interaction included, reproduces the relative energy of £1 = | and ft — | Nilsson states, in 
the N=3 major shell, and the experimental value for the first collective 2 + state. The d 
and k parameters for 154 Sm and its double beta partner 154 Gd were taken equal to those 
used in a previous publication |l3], to describe the Ml states of the mother nucleus. As 
for the remaining nuclei considered in this paper, the corresponding d and k parameters are 



the same as in Ref.|2a |26J] where one of us (A. A. R. ) described phenomenologically the 
spectroscopic properties of the major rotational bands. 

The BCS calculations have been performed within a single particle space restricted so that 
at least the states from the proton and neutron major open shells are included. Although 
the single particle energies depend on deformations here we keep calling major shell a set of 
states having the same quantum number n, according to Eq.(2.7). This truncation criterion 
defines an energy interval for single particle states. Of course, due to the level crossing 
caused by deformation, also states from the lower proton and upper neutron major shells, 
lying in the energy interval defined before, are included in the single particle space. Since 
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only the proton-neutron pair of states characterized by the same orbital angular momenta, 
participate in single beta decay processes, the single particle spaces for proton and neutron 
systems are to be the same. It is well understood that the corresponding energies for 
protons and neutrons are however different from each others for heavy isotopes, due to the 
mass dependence of the single particle mean field strength parameters (4.1). 

Pairing strengths have been fixed so that the mass differences of the neighboring even-even 
nuclei are reproduced. The results are listed in Table 1. Their values may be interpolated 
by a linear function of 1/A both for the mother : 

Gp = 1^ + 0.06931, G n = ^^ + 0.11266, (4.2) 

and the daughter nuclei: 

Gp = 1^ + 0.06765, G n = Hp? + 0.13455. (4.3) 

Slight deviations from these rules are registered for G p of 76 Ge and 76 Se and G n of 76 Ge 
and 238 Pu. It is interesting to note that although the single particle basis is different from 
the ones currently used in the literature, the results for the interaction strength is quite 
close to the standard ones. For example for 150 Sm the above equations are equivalent to 
G p = 22.86/A, G n = 25.624/ A. 

As for the proton- neutron two body interactions, their strengths were taken as in Ref . [ll | 
although the single particle basis used therein, is different from ours: 

X = ^MeV, Xl = ^MeV. (4.4) 

The A dependence for the ph interaction strength has been derived by fitting the position 
of the GT resonance for 40 Ca, 90 Zr and 208 Pb. The strength xi has been fixed so that the 
beta decay half lives of the nuclei with Z < 40 are reproduced. Certainly, the A dependence 
for x an d Xi depends on the mass region to which the considered nucleus belongs as well 
as on the single particle space. As a matter of fact our results for 76 Ge and 82 Se show that 
larger values for the ph interaction strength improve the agreement with experimental data. 
Moreover, our comparison suggests that a certain caution should be taken when the mass 
dependence given by Eq. ()4.4j) is considered as in these nucleibe x an d Xi parameters cannot 
be fixed by single beta decay half lives as it is usually done. Once the parameters defining 
the model Hamiltonian are fixed, the pnQRPA equations have been solved for the mother 
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and daughter nuclei and the output results have been used in connection with Eq. (j3.15j) 
to calculate the GT amplitude of the 2z//3/9 process. In the next step, the equation (j3.14|) 
provides the double beta half life. The phase space factor F does not depend on the structure 
of the nuclear states and therefore we take it as given in Ref.Q,!^?]]. The values for F, used 
in the present paper, correspond to g& = 1.254. Results for Mqt and half lives (T1/2) are 
given in Table 2. There, we also give the strength of the ph and pp interactions produced 
by Eq. (j4.4|) . For comparison we present also the available experimental data 



as well as the results of Ref. 



22, 



31 



32 



28| for Ti /2 . 



Nucleus 


d 


k 


G p [MeV] 


G n [MeV] 


X [MeV] 


gpp 


76 Ge 
76 Se 


1.9 


7.1 


0.300 
0.295 


0.295 
0.285 


0.35 


0.112 


82 Se 
82 Kr 


1.6 


3.5 


0.150 
0.210 


0.160 
0.215 


0.35 


0.112 


148 Nd 

148 Sm 


1.555 


10.81 


0.118 
0.120 


0.200 
0.220 


0.15733 


0.11154 


i50 Nd 
150 Sm 


1.952 


9.89 


0.160 
0.190 


0.150 
0.190 


0.15586 


0.11154 


154 Sm 
154 Gd 


2.29 


5.58 


0.190 
0.145 


0.134 
0.138 


0.15302 


0.11154 


160 Gd 

160 Dy 


2.714 


4.384 


0.160 
0.155 


0.155 
0.160 


0.14898 


0.11154 


232 Xh 
232 


2.51 


4.427 


0.120 
0.090 


0.183 
0.225 


0.11486 


0.11154 


238 |j 
238p u 


2.62 


4.224 


0.130 
0.165 


0.165 
0.235 


0.080 


0.11154 



TABLE I: The pairing and Gamow Teller interactions strength are given in units of MeV. The 
ratio of the two dipole interaction ( particle-hole and particle-particle) strengths, denoted by g pp 
is given. The list of the deformation parameter d and the factor k of the transformation (2.2) are 
also presented. The manner in which these parameters were fixed is explained in the text. 



Before discussing in extenso the results from Table 2, it is instructive to show the results 
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Nucleus 


X 


gpp 


M GT 


T 1/2 N 


present 


exp. 


Ref.[23j 


76 Ge ^76 Se 


0.35 
0.35 
0.25 


0.112 
0.112 
0.11154 


0.222 
0.149 c ) 
0.270 


5.9 • 10 20 
1.32 • 10 21 
4.05 • 10 20 


9.2+°;I • 10 20 a ) 

L1 +0.6 . 1Q 21 b) 


2.61 • 10 20 


82 ge ^82 Kr 


0.35 
0.16 


0.112 
0.108 


0.096 
0.135 


0.963 • 10 20 
0.49 • 10 20 


1.1+0-8 . 1Q 20 d) 

1.0 ± 0.4 • 10 20 e ) 
1.3 ±0.05 • 10 20 f ) 


0.848 • 10 20 


148 Nd ^148g m 


0.157 


0.112 


0.392 


2.327 • 10 19 




1.19 • 10 21 


150 Nd ^150g m 


0.156 
0.156 


0.11154 
1.50 


0.350 
0.040 


2.630 • 10 17 
1.98 • 10 19 


> 1.8 • 10 19 g ) 


1.66 • 10 19 


154 Sm -^ 154 Gd 


0.153 


0.11154 


0.327 


8.760 • 10 20 




1.49 • 10 22 


160 Gd ^160 Dy 


0.149 


0.11154 


0.170 


2.013 • 10 20 




2.81 • 10 21 


232^ _^232tj 


0.11486 


0.11154 


0.123 


4.240 • 10 21 




4.03 • 10 21 


238-y ^238p u 


0.080 
0.11282 
0.08 


0.112 
0.11154 
0.0 


0.166 
0.139 
0.171 


2.375 • 10 21 
3.340 • 10 21 
2.249 • 10 21 


(2.0 ±0.6) • 10 21 h ) 


0.914 • 10 21 



TABLE II: The Gamow- Teller amplitude for 2uf3f3 decay, in units of MeV 1 and the corresponding 
half life (Ti/ 2 ) are listed for sever al g round to g round transitions. The experimental half lives for 
the transitions 76 Ge^ 76 Se (^Ref.Q^Ref. 29] ), 82 Se^ 82 Kr ( d ) Ref. 33], e ) Ref. 32], f ) Ref.Q), 
iso Nd ^i50 gm (- g ) Ref and_^ 38 U ^ 238 Pu ( h ) Refs. (21 , 22, 3, 32|) are also given. In the last 



column the results from Ref. 28], are given. The parameters x an d Qpp are also given, 
these two cases the mother and daughter nuclei have different deformations, namely di = 1.6 and 
df = 1.9. The parameters x an d g PP are listed in the second and third columns. 



concerning the single beta decay properties of the mother and daughter nuclei. Thus, in 
Figs. 2-9 the beta minus strength of the mother and the beta plus streghts of the daughter 
nuclei, folded with a gaussian having the width equal to 1 MeV, are plotted as function of 
energy. For the lightest two nuclei, the experimental data are also presented 
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For pedagogical reasons, for these two nuclei two different ph interaction strengths are 
alternatively used. In this way one clearly sees that increasing the strength x, the transition 
strength is moved to the higher state. Thus, although the peak positions remain the same, 
since the BCS data are not changed, the first peak looses height while the second one is 
augmented. The agreement with experimental data is reasonably good. In both mother 
nuclei the center of the GT resonance is a little bit shifted, backward for 82 Se and forward 
for 76 Ge. Due to the fragmentation effect caused by the nuclear deformation, the theoretical 
result for the GT resonance has a shorter, otherwise broader peak than the experimental 
one. In order to see that a broad peak in the folded beta strength plot means, indeed, 
a fragmentation of the strength distributed among several pnQRPA states, we show the 
unfolded strength for 82 Se and 232 Th isotopes in Figs 10 and 11, respectively. For example, 
in the case of 82 Se, the folded strength exhibits a first fat peak which has a very short 
maximum before and a "shoulder" on the descending part (see the middle panels of Figs. 3). 
From Fig. 10 one sees that to these details correspond pnQRPA states which carry a strength 
represented by sticks which dominate the grass spread around. 

The total strength of the GT resonance is about the same as the corresponding exper- 
imental data. However the fragmentation is causing a broad resonance which results in 
having a shorter peak. This fact might raise the question whether the deformation consid- 
ered in the present paper is too large. Indeed, the neutron system of 82 Se is almost spherical 
since N (the neutron number) is close to a magic value. Due to this feature we repeated the 
calculations with a very small g? (=0.2) which is close to the spherical limit. As seen from 
the right panel the height of GT resonance corresponding to the new deformation is close to 
the experimental result. Since we kept the same parameter for the single particle mean field 
parameters, e.g for the parameter k defined in Eq. (2.2), the theoretical curve is shifted by 
about 1 MeV with respect to the experimental data. Of course, the position of the GT res- 
onance depends also on the pairing strengths. The parameters mentioned above where kept 
the same as for the initial deformation case (d=1.9), in order to judge, by comparison, the 
effect coming from the nuclear deformation. In conclusion, going from deformed to spherical 
single particle basis the GT resonance peak is getting higher and the width narrower. Actu- 
ally when the calculated strength distribution is compared with the experimental data one 
has to restrict the discussion only to the position of the GT centroid and the total strength, 
since there is no experimental information about the resonance width. The narrow whith 
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seen, however, in Figs. 2 and 3 is caused by folding a single number indicating the total f3~ 
strength for the GT resonance which has the centroid at a given energy, with a Gaussian 
having the width equal to 1 MeV. 

From the folded beta minus strengths graphs, we see that for heavier double beta emitters 
there exists a small peak lying above the GT resonance. 



Nucleus 


1st peak 


2nd peak 


3rd peak 


Transition 


Strength 


Transition 


Strength 


Transition 


Strength 


76 Ge 


*(*f§3) 


1.084 


K44i)^7T(4d§§) 

iyC3 f lh) -> 7r(3fM) 

V-V 2 2/ v ^ 22/ 


1.740 
3.371 


K3/5§) -tt(3/§§) 


1.718 


82 Se 


K3/§§)-7T(3/5§) 


1.554 


i/(3p§§)-7r(3p§±) 


4.501 


i/(3p§i-7r(3p§§) 
^(3/11)^71(3/11) 


1.686 
3.310 


148 Nd 


"\^y 22) " \ y 2 2 J 


1.077 
1.141 


K4 5 ||)^7r(44|) 


1.453 
12.437 
1.015 






150 Nd 


\ y 22/ v y 22/ 


1.901 


V" 2 2/ v 2 2 / 

K44|)-vr(44f) 
K44|)^7r(4 5 |f) 


1.246 
1.178 
1.087 
7.647 


i/(4<Z§i1 -> 7r(4eZ|§) 

V "2 2/ V 2 2/ 

K44|)-7r(44f) 


5.386 
2.370 


154 Sm 


i/(5/5§) -»7r(5/5§) 

I/(6if |)-»7T(6if |) 

K4<7l|)-vr(44i) 


1.220 
1.031 
3.362 


i/(4(i|i) -► 7r(4d§§) 


4.214 
2.088 


K44|)-7t(44|) 

^(441)^^(441) 


1.537 
4.380 


160 Gd 


^fD-^iif) 

K4 5 H)-vr(4 5 ||) 


2.222 
1.622 


^(441)^77(441) 
i/(6zf |) -► tt(6^|) 


2.551 
1.362 
9.922 


^(441)^^(441) 
^441)^^(441) 


1.028 
1.204 



TABLE III: Here are listed the strengths carried by the pnQRPA states contributing to the first, 
second and third (if any) peaks from the upper panels of Figs. 2-9. On the left side of these numbers 
the 2qp configurations closest in energy to the corresponding pnQRPA states, are given. Actually 
this is the dominant configuration of the chosen pn phonon state. 
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Nucleus 



1st peak 



Transition 



Strength 



2nd peak 



Transition 



Strength 



3rd peak 



Transition 



Strength 



232 



v(5h 



11) 

2 2> 



tt(5/ 



5 3\ 
2 2/ 



1.810 
1.569 





— * 


7T(4d|i) 






/r 7 11 5 A 


/T7 11 3\ 




fri 9 5\ 

7r(5n||) 


/r r 11 5 \ 

i/(5/i^§) 




/r 7 9 3 \ 
7r(5/l|§) 






c/i ■ 11 7\ 


/ r , 11 7\ 




( f r 11 5 \ 

t(5/t§) 


u(6if |) 




//-. • 13 1 A 


■ 13 3 \ 




(n ■ 11 5 \ 


i/(5/i§|) 




vr(5/i§±) 


i/(5/i§±) 




tt(5/i§§) 






tt(5/i§|) 






tt(44§) 












*(5/§§) 



1.900 
9.913 
7.157 
2.190 
2.392 
8.709 



i/(5p§±) 



ir(5p 



3 3\ 
2 2/ 



2.548 



238 



u 



51) 

2 2/ 



*(5/§§) 



2.253 
1.079 
li 



4.259 
2.868 
1.112 
6.201 
2.912 
2.019 
3.662 
1.222 



i/(5pf±) 



7r(5p§±) 



1.397 



TABLE IV: Continuation of Table III. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 




pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


76 Ge 


7.033 


1.084 


10.850 


1.740 


12.602 


1.718 








11.605 


3.371 






82 Se 


6.939 


1.554 


10.920 


4.501 


11.701 


1.686 












12.291 


3.310 


148 Nd 


9.397 


1.077 


12.028 


1.453 








10.047 


1.141 


12.269 


12.437 












12.429 


1.015 






i50 Nd 


8.600 


1.901 


11.263 


1.246 


12.939 


5.386 








11.531 


1.178 


13.217 


2.370 








12.281 


1.087 












12.597 


7.647 






154 Sm 


10.475 


1.220 


11.986 


4.214 


13.189 


1.537 




11.047 


1.031 


12.696 


2.088 


13.303 


4.380 




11.434 


3.362 










160 Gd 


10.748 


2.222 


12.457 


2.551 


15.334 


1.028 




11.163 


1.622 


12.857 


1.362 


15.850 


1.204 








13.369 


9.927 







TABLE V: The energies of the pnQRPA states which give the largest strength contributions to 
the peaks in Figs. 2-9, upper panels. The carried strengths are also given. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 




pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


232 rp^ 


12.895 


1.810 


15.622 


1.900 


19.559 


2.548 




13.578 


1.569 


15.952 


9.913 












16.367 


7.157 












16.593 


2.190 












16.731 


2.392 












16.942 


8.709 






238 |j 


13.641 


2.253 


16.079 


4.259 


21.137 


1.397 




14.792 


1.079 


16.306 


2.868 








15.219 


1.888 


16.452 


1.112 












16.559 


6.201 












16.613 


2.912 












16.831 


2.019 












17.391 


3.662 












18.232 


1.222 







TABLE VI: Continuation of Table V. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 


Transition 


Strength 


Transition 


Strength 


Transition 


Strength 


76 Se 


7T(3/|5)- i/(3/§§) 


0.313 


^(3/75)^(3^5) 
^(3/75)^(3^3) 

7r(3/||) ^(3/||) 


0.116 
0.027 
0.021 


tt(3/§§)- i/(3/5§) 


0.013 


82 Kr 


vr(3/|l)^K3/||) 


0.135 


^(3 p 33)^^( 3p ll) 


0.118 


^(3/73)^^5 5) 


0.019 


148 Sm 


tt(5^|) ^v{hh%\) 

VT^/l^) _^ ^(5^93) 
V 2 2' V 2 2 / 

tt(5^|) i/(5h§|) 


0.139 
0.151 
0.165 


vr(44i)^K44l) 

7r(4o||) ^ j,(4 g 77) 


0.037 
0.017 


^(441)^^(441) 

7r(4o|§) _^ ^(4095) 
V i> 1 1 ) \H HI 


0.011 
0.012 


150 Sm 


tt(5^|) -► i/(5h§|) 
^(5^1) i/(5/i§|) 
7T(5^|) ^u{bh%\) 


0.227 
0.134 
0.103 


vr(5/i|i) i/(5/i§§) 

V 22^ \ 2 2 / 


0.012 
0.018 






154 Gd 


7r(4d||)^K4d||) 
^(5^1)^^(541) 


0.081 
0.140 


^113)^^93) 

^(5^1)^^(541) 


0.061 
0.046 


^113)^(5^9 1) 

7T(4d§I)^I/(4d||) 

^(409 7) ^ ;,(4n7 7) 

V i» 2 2 / V j 2 2 / 


0.039 
0.019 
0.016 


160 Dy 


7r(4d§§) -»• i/(4d||) 
7r(5^|)^K54|) 

2 2/ V w 22/ 


0.050 
0.076 
0.365 


7T(5/l§i) -»• ^(54f) 
^99)^^77) 

7r(6i¥f) _^ !/(6iiJi) 

\ 2 2/ ^ 2 2/ 

tt(5^§) -»• ^(54f) 


0.067 
0.113 
0.059 
0.073 


rr(6zf i) i/(6i^§) 

7T(4d§l)^I/(4d||) 

\ 2 2 / V 2 2 / 


0.063 
0.051 
0.054 


232 


vr(6^|) i/(6i^§) 

V 2 2' v 11' 
7T(6zf §)->l/(6i^ |) 

,r(6i^§)-i/(6i^§) 


0.011 
0.013 
0.011 


tt(6^§) i/(6i^§) 

V 2 2' V 2 2' 


0.031 


vr(6if |) i/(6i4±§) 

\ 2 2' V 2 2' 


0.008 


238p u 


^117)^^(6^7) 
vr(5/||)^K5/||) 


0.008 
0.009 




0.010 


7r(6^)^(6if!) 


0.005 



TABLE VII: The same as in Table III but for the lower panels of Figs. 2-9. 
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Nucleus 


1st peak 


2nd peak 


3rd peak 




pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


pnQRPA energy 


Strength 


76 Se 


2.684 


0.313 


4.904 


0.116 


13.028 


0.013 








6.169 


0.027 












6.627 


0.021 






82 Kr 


2.158 


0.135 


5.313 


0.118 


8.172 


0.019 


148 Sm 


2.196 


0.139 


6.209 


0.037 


10.029 


0.011 




2.301 


0.151 


7.408 


0.017 


11.899 


0.012 




2.427 


0.165 










150 Sm 


3.051 


0.227 


6.420 


0.012 








3.102 


0.134 


7.369 


0.018 








3.419 


0.103 










154 Gd 


2.071 


0.081 


3.169 


0.061 


3.928 


0.039 




2.491 


0.140 


3.600 


0.046 


4.170 


0.019 












4.401 


0.016 


160 Dy 


4.587 


0.05 


5.616 


0.067 


6.209 


0.063 




4.756 


0.076 


5.742 


0.113 


7.469 


0.051 




5.392 


0.365 


5.830 


0.059 


7.559 


0.054 








5.968 


0.072 






232 


3.107 


0.011 


9.531 


0.031 


16.069 


0.008 




4.917 


0.013 












6.656 


0.011 










238p u 


3.378 


0.008 


11.751 


0.010 


17.689 


0.005 




3.415 


0.009 











TABLE VIII: The same as in Table V, but for the lower panels of Figs. 2-9. 
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Of course, this feature is mainly caused by the fact that while for Ge and Se transitions 
the states contributing to the GT resonance have an energy separated by a gap from the 
upper two qp dipole configurations, in the heavier nuclei such energy gap does not exist due 
to both deformation effect on single particle energies and the fact that the last filled state 
is far away from a major shell closure. 

In Tables III and IV we list the single particle f3~ transitions characterized by the fact that 
the corresponding two quasiparticle energy is the closest one to a pnQRPA state contributing 
to the n-th peak with the strength given at its right side. The pnQRPA energies for the states 
bringing the strength listed in Tables III and IV are given in Tables V and VI, respectively. 

The single particle transitions v(nljl) — > 7r(nlj'I') which coherently contribute to the 
collective transition + — > 1 + , are characterized by the change of quantum numbers j and I 
by at most one unit i.e. |Aj|=0,l and |A/|=0,1. The (Aj, AI) values for the single particle 
transition which represents the dominant component of the pnQRPA state which carries 
the maximal strength in a GT resonance are (1,1) ( 76 Ge), (0,l)( 82 Se), (0,1) ( 148 Nd), (0,1) 
( 150 Nd), (1,0) ( 154 Sm), (1,1) ( 160 Gd), (1,1) ( 232 Th), (0,1)( 238 U). From Tables III and IV, it 
results that the GT resonances are admixtures of AJ + Aj = 1,2 transitions. 

A common feature for all nuclei considered in the present paper is that the dominant 
component of the pnQRPA state, i.e. that component which is excited with the largest 
probability by the GT transition operator, involves single particle states with small I. In- 
deed, such transitions vl — > ixV have either I or /' equal to \ or |. In the cases of 154 Sm the 
angular momenta, involved in the transition are equal to each other. The common value is |. 
In Figs. 2-9 we give also the folded 2qp strength. One notes that the pnQRPA correlations 
push the strength to the higher energy. One of the main effects, comparing it with the 2qp 
picture, is that it concentrates most of the strength in the GT resonance which is the most 
collective pn excitation in the intermediate odd-odd nucleus. 

Now, let us focus our attention on the f3 + strength distribution in the daughter nuclei. 
These strengths are much weaker in magnitude than those characterizing the (3~ strength 
in the mother nuclei. Another difference between the two processes is that for the single 
beta minus process the maximum strength is concentrated at relatively high energy, around 
the Gamow- Teller resonance, while in the beta plus decay, most of the strength lies around 
5 MeV. Indeed, from the lower panels of Figs. 2-9, one sees that the first peak, is the 
highest one. Exception from this rule is 232 Th where the (3 + strength is quite small and 
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its distribution has a peak lying around the energy of 10 MeV. Switching on the QRPA 
correlations one notices a decrease of the 2qp strength. Actually the difference in strengths 
which appear for the peaks is distributed among the remaining pnQRPA states, the total 
amount of strength in the two pictures being the same. The two qp configurations which 
contribute most to the first peaks shown in Figs. 2-9, are listed in Table VII. They are the 
dominant components of pnQRPA states with the energies given in Table VIII. As in the 
case of the single f3~ decay of the mother nuclei, here also most of the dominant transitions 
take place between states of low angular momenta (/ = |, ~). However, due to the small 
magnitude of the transition strengths here one notices several transitions between states 
with angular momenta equal to |, |, §. 

Inspecting the expression of the double beta transition amplitude, one notes that the 
numerator of a chosen term from the sum, has three factors: i) one which determines the 
strength of the f3~ transition to a particular state ii) one whose hermitian conjugate 

matrix element describes the (3 + transition to a state \lk')f and iii) the overlap of the 
states reached by the decay of the initial and final nuclei, respectively. One expects that 
the maximum contribution to the GT transition amplitude Mqt is achieved when the two 
single beta matrix elements are maximal and moreover the overlaps of the dipole states in 
the odd-odd system is maximum. Therefore, we could ask ourselves whether among the 
peaks in the upper panels and those of the lower panels there are pairs of peaks determined 
by states of maximal overlap. From Tables III and VII one could identify many such pairs 
of peaks from the beta minus and beta plus strength distributions. For illustration we 
mention only one example. In 148 Nd the maximum contribution to the GT resonance is 
brought by the pnQRPA state of energy equal to 12.269MeV. Indeed, the corresponding 
strength is 12.437 and moreover this is the leading strength. The dominant amplitude for 
this state corresponds to the single particle transition i/(4g||) — > 7r(4g||). On the other 
hand the strength distribution for 148 Sm shows a third peak determined by the pnQRPA 
state of energy equal to 11.899 MeV. Since the dominant two quasiparticle component of 
this collective state, corresponds to the single particle transition 7r(4g||) ~^ K^ff)' as 
shown in Table VII, one expects that the overlap of the pnQRPA states mentioned above is 
maximally large. 

The single beta decays strengt 
under the name of Ikeda sum rule 
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is of a given nucleus satisfy the N-Z sum rule, known 
. Our predictions for j3~ and j3 + strengths satisfy the 



Ikeda sum rule in the heavy isotopes while for 76 Ge and 82 Se small deviations of 3 and 1.7% 
respectively, are registered. 

Let us analyze now the results for the double beta process, given in Table II. The dipole- 
dipole interaction strengths have been chosen as given by the empirical formula (J4.4j) . The 
ratio of the pp and ph interaction strengths determines the g pp factor. As we already men- 
tioned before, this A dependence for the interaction strengths depends on single particle 
basis as well as on the truncation of the single particle space and therefore its validity for 
the present formalism is questionable. As a matter of fact for the lightest nuclei a larger 
value for x approaches better the experimental situation while for 238 U a smaller value is 
more suitable. Definitely, the safe way of fixing the ph and pp interaction strengths would 
be to fit the position of GT resonance centroid of the odd-odd intermediate nucleus and the 
half-lives of the (3 + decay of the of unstable nuclei in this mass region, respectively. However 
since for the cases considered here the experimental data mentioned above are lacking , we 
adopted the empirical formula ()4.4|) just to obtain some reference results to be compared 
with the ones obtained with the same interaction but different single particle basis. 

The denominator from the equation defining Mqt involves the Qpp values and the ex- 
perimental energy of the first 1 + . These values are given in Table IX. Except for the case 



Nucleus 


76 Ge 


82 Se 


148 Nd 


150 Nd 


154g m 


160 Gd 


232 


238 |j 


Q^[MeV] 


2.04 


3.01 


1.93 


3.37 


1.25 


1.73 


0.85 


1.15 


Nucleus 


76 As 


82 Br 


I48p m 


150p m 


154 Eu 


160 Tb 


232 Pa 


238 Np 


E 1+ [keV] 


44 


75 


137 


137* 


72 


139 


1000* 


244 



TABLE IX: The Q^/3-values for mother nuclei are given in units of MeV. In the lowest row, the 
experimental energy for the first 1 + states in the intermediate nuclei are given in units of keV.*) For 



150 p m an( j 232 p a thgj-g are no t available data. Therefore we take the ad hoc values characterized 
by an asterisk. Actually changing E 1 + within 1 MeV does not modify the order of magnitude of 
T1/2 

of ( 76 Ge; 76 Se) all other pairs of (mother;daughter) nuclei are characterized by only slightly 
different nuclear deformations. For this reason in our calculations the nuclear deformations 
of mother and daughter nuclei have been considered equal to each other. The results for 
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Mqt and T1/2 are listed in Table II. They are compared with the experimental available 



data as well as with the predictions of those of Ref. 



28j. Table II shows a good agreement 



between the predicted T-i /? for 82 Se and 238 U, and the corresponding experimental half life 
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3,0 



and Ref. 



21 



22, 



331 ] . respectively. Our prediction for the half life of 



given in Ref. 

150 Nd is 69 times lower than the corresponding lower experimental limit. As shown in the 
second row for this nucleus, this discrepancy can be recovered by changing g pp to a value 
equal to 1.50. In this context we would like to mention that while for the lightest two nuclei 
from Table II, the Mqt function of g pp shows a very abrupt decreasing part, for the heavier 
nuclei the cancellation point is reached with a curve of a moderate slope. In the case of 
150 Nd the cancellation value of g pp is larger than 1.8 and therefore the adjusted value of 1.5 
is still far away from the critical point where the pnQRPA breaks down. 

The predicted half life of 76 Ge shown in the first row of Table II is only slightly smaller 
than the lower limit of the corresponding experimental data. For this case, however, the 
deformations for mother and daughter given in Ref. [40] are quite different from each other. 
This feature challenged us to consider in our calculations different deformations for 76 Ge and 
76 Se. Therefore, we repeated the calculations for the decay of 76 Ge with the deformation 
d m = 1.6 for the mother and da — 1.9 for the daughter nucleus. The pairing strengths for 
the new value for the nuclear deformation acquired by the mother nucleus are G p =0.290, 
G n =0.280. The results are shown in the second row of Table II. From there one remarks a 
good agreement with the experimental data for T\ji- 

To conclude, considering different deformations for mother and daughter nuclei decreases 
the overlap matrix elements involved in Mqt- Due to this effect the T1/2 value is increased. 
It is an open question whether considering different deformations for 150 Nd and 150 Sm would 
wash out the big discrepancy with the experimental data shown in the first row for the decay 
of 150 Nd. At the first glance one may say that in order to have a positive answer to the 
question formulated above, one needs a larger difference between the two deformations than 
indicated in Refs. [39,40]. Since the half life is sensitive to the pairing properties one may 
suspect that for this case the proton-neutron pairing might play an important role. 

Comparing the results for Tu 2 obtained with our formalism with those obtained in Ref. 
by a different method one notices that for four emitters, 148 Nd, 150 Nd, 154 Sm, 160 Gd, our 
predicted half lives are shorter than in the above quoted reference while for the remaining 
nuclei the ordering of the half lives is opposite. In some cases the difference between the two 
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sets of predictions are in the range of two orders of magnitudes. Since the two methods are 
based on different approaches for the transition amplitudes, it is an open question whether 
these big discrepancies could get a consistent justification. 

Finally we addressed the question of how the GT transition amplitude depends on g pp 
and whether this dependence is influenced by the nuclear deformations. The results of our 
investigation are presented in Figs. 12, 13. The input parameters of single particle states 
and pairing interactions are those of 76 Ge. From Fig. 12 one sees that for small values of 
9pp(< 0.5), M-gt depends monotonically on d while for g pp > 0.5 this property is lost. The 
repulsive character of the pp interaction causes the cancellation of Mqt for a g pp around 3. 
As seen from Fig. 12 the cancellation point depends on deformation. Also the curves look 
of Fig 12 is not changed, the value of g pp where M GT is canceled is quenched by the factor 
by which the strength \ is increased when one passes from Fig. 12 to Fig.13. In this context 
we recall that for spherical nuclei, the cancellation, corresponding to the value of x which 
reproduces the position of the GT resonance, takes place for g pp ps 1. For this value the 
relation between the matrix elements of the ph and pp two body interactions is given by the 
Pandya transformation. From Fig. 13 it results that the cancellation points depend, as we 
already said, on deformation. It is an open question whether the deformation dependence 
of the GT resonance is such that the cancellation point of Mqt is always brought to about 
1. 



V. SUMMARY AND CONCLUSIONS 



The main results described in the previous sections can be summarized as follows. The 
two neutrino double beta decay transition amplitudes and half lives for eight isotopes have 
been calculated within a pnQRPA approach based on a projected spherical single particle 
basis. 

The single particle energies are approximated by averaging a particle-core Hamiltonian on 
the projected basis. Due to the fact that the core volume conservation is properly taken into 
account, the resulting energies depend on deformation in a similar manner as Nilsson levels. 
This feature suggests that the results for two neutrino double beta decay rate provided by 
a pnQRPA formalism with such a basis will be essentially different than those obtained in 
Ref.|9j where single particle energies depend linearly on deformation. 
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First we studied the (3~ and (3 + strength distributions for mother and daughter nuclei 
respectively. Both types of strengths are fragmented due to the nuclear deformation. The 
position as well as the width of the GT resonance depend on nuclear mass. Moreover, while 
the GT resonance lies in the upper part of the pnQRPA energy spectrum (the meaning 
of this statement is that beyond the GT resonance there is only a little strength left to 
be distributed) the highest peak in the folded strength distribution plot for the f3 + decay 
appears always for low energies. This feature suggests that the GT resonance is mainly 
influenced by the ph while the peak in the (5 + strength distribution, by the pp channels 
of the dipole-dipole two-body interaction. It seems that there is a correspondence between 
the pnQRPA states of mother and daughter nuclei contributing most to the folded strength 
distributions. The associated states, due to the correspondence mentioned above, have 
maximal overlap and therefore give the main contribution to the Mqt value. From Figs. 
2-9 one remarks that the GT resonance strength depend on the atomic mass. The larger 
A, the larger the height of the resonance. Moreover the energy of the resonance center is 
also an increasing function of A. For two emitters, 82 Se and 154 Sm, the GT resonance has 
a doublet structure. This reminds us of the doublet structure of the dipole giant charge 
preserving resonances due to the coupling to the quadrupole degrees of freedom. Actually 
in this case also the doublet structure is a deformation effect and by this an effect caused 
by the quadrupole coordinates of the core. 

The Mgt and Ti/ 2 values have been first calculated by considering equal deformations 
for mother and daughter nuclei. The A dependence of the ph and pp proton-neutron dipole 
interaction is taken as in Ref . \v\ . The agreement with experimental data concerning the 
Tx/2 value of 82 Se and 238 U is very good good. The result for 76 Ge is slightly smaller than 
the experimental data. The discrepancy was removed by considering the deformation for 
76 Ge different from that of 76 Se. Indeed, this is the only case where according to Refs. [30,40] 
the deformations for the two nuclei involved in the double beta decay are quite different. 
To bring the theoretical value of Ti/ 2 for 150 Nd in agreement with the experimental data 
one needs a larger deformation difference than given in literature [39,40]. Moreover, the 
pairing strength should deviate very much from what the difference of neighboring even-even 
isotopes masses requires. Due to this feature for this case we reproduced the experimental 
half life by changing g pp from 0.11154 to 1.5. Note that the critical value of g pp for this 
isotope is 1.8. It is noteworthy that for isotopes, 76 Ge, 82 Se, 238 U, where the calculated 
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half lives agree with the corresponding experimental data the values used for g pp are small 
which results in having a small effect coming from the pp interaction on this observable. As 
shown in Table II for 238 U, cancelling g pp does not alter the agreement with the experimental 
data. Then it arises the question whether the pp interaction is really needed at all in order 
to describe quantitatively the double beta decay process. Is the large sensitivity of the 
single (3 + matrix elements, pointed out for the first time by Cha in Ref. 4l|, a real effect 
or just an artifact caused by the instability of the pnQRPA ground state 8|? As shown in 
Table II for 76 Ge, taking different deformations for mother and daughter nuclei brings an 
important effect on T1/2 but not a dramatic change as claimed in Ref. [42]. The difference 
between the two descriptions consists in the fact that here the overlap matrix elements of 
the states in the mother and daughter nuclei are estimated in a manner consistent with the 
pnQRPA approach, while in the quoted reference the phonon operators are dissociated and 
the overlaps are calculated within the BCS and particle representations. Of course, in the 
latter case it is not possible to get a real hierarchy of the effects involved. 

As we stressed in Ref . jiol] . going beyond pnQRPA approach, some forbidden processes 
might become possible. As an example, we studied the double beta transition on excited 
collective states. In the near future we shall investigate whether increasing the deformation 
for the daughter nuclei would substantially increase the reduced decay probability for such 
processes. 

VI. APPENDIX A 

In order to calculate the overlap matrix which enter the Mqt expression, we have to 
express the phonon operator for the mother nucleus in terms of the phonon operator of the 
daughter nucleus, following the boson expansion procedure 

rU<, k) = k')rl(f, k 1 ) + z{(k, fcOiW/, AOC-) 1 """]. (A-i) 

k< 

where the amplitudes W and Z can be easily calculated as follows: 

w/(k,k') = t <o| [r^c/, AO , rl^Ci, *)] [o>/, 

Z/(M0 = l (0\[Tl(f,k)A,-^k')(-) 1 -»]\0) f . (A.2) 
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It is clear that once these amplitudes are calculated, the overlap matrix elements are readily 
obtained: 

i <l*|l , fc >/ = W/(A; > * , )<<0|0> / (A.3) 

provided that the overlap of the two vacua is known. In what follows we shall describe the 
necessary steps to derive the expressions of the two factors from the right hand side of the 
above expression. By a direct calculation one finds: 



W/(k, k')* = (X k (i;pn)X k/ (f;p'n') ~Y k ( i -,pn)Y k if-,p'n%(0\[A 1 ^ ] pn),Al l (f- 1 p'n')]\0) f 
+ (X k (i;pn)Y k if;p'n')-Y k (t;pn)X k if;p'n%(0\l(^ 

(A.4) 

The symbol "*" stands for the complex conjugation operation. The matrix elements of 
the commutators of the two quasiparticle operators are expressed further in terms of the 
anti-commutator of the single particle operators which is calculated as explained in Ref . : 

{c a iM^ci PM ,(f)} = AcyoAc;,^ x 

j 

x ( j c) (i, f)5 ir 5 jf 5 MM > = 0^5 n t5 j f5 aa >5 M M' ( A -5) 
Here Nj 9 ^ denotes the norm of the core projected state: 

A{d) = N { f\d)Pi I0 e d ^o-^) lo)b . . (A.6) 

where |0)& denotes the vacuum state for the quadrupole bosons. The overlap matrix for the 
initial (i) and final (f) core states is denoted by Of (i, /), and has the expression: 



of{iJ) = (<p JM (d)\<pj M (d?)) ^ N?(d)N^(d')e-^ L (2L + l)lf\Vddf). (A.7) 

where the factor I/is defined by Eq. ()2.1ip . The initial nucleus deformation is denoted by 

d while d! stands for the deformation parameter characterizing the final nucleus. Note that 
in the limit d' — *> d we have 

0f>(i,f)^l, , (A.8) 
With these details, one further obtains: 
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(A.9) 

and 



i(0|[(-) 1 -M t li _ /1 (i;pn),Al /1 ,(/;p'n')]|0> / = ^^^(^^-^^(EC^-^^O^O. -''' 



(A.10) 



Thus, the final expressions for the amplitudes W and Z are: 



k') = Yl [(X k (i;pn)X k ,(f;pn) - Y k (i;pn)Y k ,(f;pn))(U;Uf + V;v/)(W n Ul + V l n Vl) 

pn 

+ (X k (i;pn)Y k ,(f;pn) - Y k (i;pn)X k if;pn))(U;V p f ~ V##)(tW - 
= ^ [(I fe (i;p)X fc ,(/;pn) - y,(i;pn)y,(/;pn))([/;y/ - ^^(^ - KJ^) (A.11 

pn 

+ (X k (i;pn)Y k ,(f;pn) - Y k (r lP n)X k if;pn))(U;uf + V;vJ)(U l n lJf + ^y/) 
As for the ground states overlap, the result is 



i <0|0> / « i(BCS\BC3) f 

= II + W (°y /p )1 II \uM + v*yt (off 1 ™)'' 



(A.12) 



If the mother and daughter nuclei are characterized by the same nuclear deformation, the 
corresponding overlap matrix elements are obtained from the above formulae by replacing 
d! by d. A good approximation of the resulting equation is given by the expression: 



i{\ k \l k i)f = Y [X k (i,pn)X k ,(f,pn) - Y k (i,pn)Y k (f,pn)} . 



(A.13) 



pn 



The drawback of the procedure described above is that both factors of Eq lA.3l are evaluated 
within the framework of the BCS approximation while the matrix elements describing the 
single j3 transitions are calculated within the pnQRPA approach. Moreover, due to the 
overlap factors Ofj, even the border of BCS frame is crossed toward the particle represen- 
tation. This inconsistency of the levels of approximations makes the method doubtful, since 
it is not possible to define an hierarchy of various effects (3). For example one could take 
care of a negligible contribution otherwise ignoring an important one. However, extending 
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the spirit of the RPA approach to the case of different deformations for the initial and final 
nuclei one obtains: 

Wf(k,kT) = ^[X fc (i;pn)X fc ,(/;pn)-y fc (i;pn)y fe ,(/;pn)], 

pn 

,(010); = 1 (A.14) 

In this way the overlap matrix and the matrix elements characterizing the initial and final 
nuclei are treated in an unitary fashion. The numerical calculations presented in the present 
paper correspond to the overlap matrix determined by Eq. (jA.14|) . We postpone for a forth- 
coming paper, the description of the (3^ matrix elements within a higher pnQRPA approach 
consistent with the procedure presented in this Appendix for calculating the overlap matrix 
elements. 
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FIG. 1: (Color on line)Proton single particle energies for the N=3 and N=4 major shells, given in 
units of huo, are plotted as function of the deformation parameter d. The quantum numbers on 
the right hand side are nljl defined in Eqs. (2.5) and (2.7) 
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FIG. 2: (Color on line)Single /3~ strength, for Ge, and single f) + strength for 76 Se, folded with a 
Gaussian function having the width of 1 MeV, are plotted as a function of the energy within the 
BCS and pnQRPA approximation, for three values of the particle-hole interaction strength, x-The 
left and middle panels correspond to d=1.9 while the right panel to = 1.6 and = 1.9. For 
X = 0.35 and x = 0-4, we also give the experimental results from Ref. 
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FIG. 3: (Color on line)The same as in Fig. 2, but for the f3 of 82 Se and the j3 + of 82 Kr. 
Experimental data are from Refs 34]. The right panels correspond to d=0.2. 
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FIG. 4: The same as in Fig. 2, but for the (3~ of 148 Nd and the (3+ of 148 Sm. 
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FIG. 5: The same as in Fig. 2, but for the (3~ of 150 Nd and the (3+ of 150 Sm. 
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FIG. 6: The same as in Fig. 2, but for the /T of 154 Sm and the (3+ of 154 Gd. 
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FIG. 7: The same as in Fig. 2, but for the 0~ of 160 Gd and the 0+ of 160 Dy. 
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FIG. 8: The same as in Fig. 2, but for the p~ of 232 Th and the 0+ of 232 U. 
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FIG. 9: The same as in Fig. 2, but for the /3~ of 238 U and the [3 + of 238 Pu. 
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FIG. 10: The single (3 , for 8 Se, and single (3 + , for 8 Kr, are plotted as a function of the energy 
within the pnQRPA approach, for two values of the particle-hole interaction strength, x- 
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FIG. 11: The single (3 , for 232 Th, and single (3 + , for 232 U, are plotted as function of the energy 
within the pnQRPA approach, for x = 0.11. 
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FIG. 12: (Color on line) The Gamow- Teller amplitude Mqt for the transition 2v(3(3 is represented 
as a function of g pp , for a particular value of the particle-hole interaction strength, x = 0.1. 
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FIG. 13: (Color on line) The same as in Fig. 12 but for \ = 0.3. 
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